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Abstract — An asymptotic technique is developed to find the 
Signal-to-Interference-plus-Noise-Ratio (SINR) and spectral effi- 
ciency of a link with N receiver antennas in wireless networks 
with non-homogeneous distributions of nodes. It is found that 
with appropriate normalization, the SINR and spectral efficiency 
converge with probability 1 to asymptotic limits as N increases. 
This technique is applied to networks with power-law node 
intensities, which includes homogeneous networks as a special 
case, to find a simple approximation for the spectral efficiency. 
It is found that for receivers in dense clusters, the SINR grows 
with N at rates higher than that of homogeneous networks 
and that constant spectral efficiencies can be maintained if the 
ratio of N to node density is constant. This result also enables 
the analysis of a new scaling regime where the distribution of 
nodes in the network flattens rather than increases uniformly. 
It is found that in many cases in this regime, N needs to 
grow approximately exponentially to maintain a constant spectral 
efficiency. In addition to strengthening previously known results 
for homogeneous networks, these results provide insight into the 
benefit of using antenna arrays in non-homogeneous wireless 
networks, for which few results are available in the literature. 

I. Introduction 

Multiple antenna systems can improve the performance of 
point-to-point wireless links by increasing data rates through 
spatial multiplexing and array gain as well improving robust- 
ness through diversity. In wireless networks, multi-antenna 
arrays can be used to mitigate interference enabling simulta- 
neous transmissions. The data rates achievable using antenna 
arrays for interference mitigation is highly dependent on 
the distribution of nodes in space which impacts inter-node 
distances, signal and interference levels, and hence data rates. 

Thus, the spectral efficiencies achievable in wireless net- 
works with spatially-distributed multi-antenna nodes have re- 
ceived much attention. Chen and Gans [1| found that in order 
to maintain constant per-link data rates in ad-hoc wireless 
networks, it is necessary to increase the number of receiver 
antennas linearly with the number of users. Govindasamy et. 
al. [|2;| found that the SINR of a representative link in a wireless 
network where signal power decays according to a power-law 
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with path-loss exponent a, grows as the ratio of the number 
of antennas N to the user density p raised to the power 
a/2. Jindal et. al. [3| found that it is possible to increase 
the area spectral efficiency of wireless networks linearly with 
the number of receiver antennas by increasing the density 
of simultaneous transmissions using a partial zero-forcing 
receiver. Hunter et. al. [4] and Kountouris and Andrews 
consider several transmit encoding approaches with the former 
considering approaches such as antenna selection, maximal 
ratio transmissions and space-time orthogonal block codes, 
and the latter analyzing dirty-paper coding, zero-forcing and 
block-diagonalization. Ali et. al. |6), and Louie et. al. Q found 
the CDF of the SINR in networks with Linear Minimum- 
Mean-Square-Error (MMSE) receivers in Rayleigh fading with 
single and multi-stream transmissions respectively. Vaze and 
Heath [8|, Hunter and Andrews |9| and Govindasamy et. al. 
[10 1 considered multi-stream transmission with channel-state- 
information (CSI) at transmitters using different assumptions 
on the transmitter and receiver processing. These works as- 
sume a uniformly random distribution of nodes on the plane 
where the intensity of the nodes does not vary with location. 
Networks where node intensities vary with location arise in 
many practical situations such as when users congregate at 
hot-spots. For non-homogeneous distribution of nodes, Ganti 
and Haenggi [ 1 1 1 develop techniques to analyze wireless 
networks with clusters of users and Tresch and Guillaud in 
[12 1 analyze interference alignment in such networks. Win et. 
al. iFPJI develop a general set of tools to analyze interference 
in wireless networks with possibly non-homogeneous node 
distributions. 

In this work we develop an asymptotic technique to find 
the spectral efficiency and SINR of a link with multiple 
antennas at the receiver, in a network with a non-homogeneous 
distribution of nodes in the interference-limited regime, with 
Linear MMSE receivers. In contrast to [11| and [12|, we 
assume that the intensity function is fixed which can arise 
in situations where the locations of hot-spots are known. The 
spatial distribution of nodes plays an important role in the 
rate at which SINR grows with the number of antennas at the 
receiver and hence strongly influences the trade-off between 




Fig. 1. Randomly distributed points with linear radial intensity 



the increased hardware costs of antenna arrays and increased 
data rates. For instance, we find that in a network where the 
intensity of users decays inversely with the distance from the 
receiver (as illustrated in Figure [TJ, the SINR grows as N a 
compared to N a / 2 for a uniform intensity of nodes. This 
result indicates that antenna arrays can provide an even higher 
increase in SINR in such networks compared to uniformly 
random networks. 

We address the added complexity of non-uniform node 
intensities by using an asymptotic analysis, taking the radius of 
a circular network and number of receiver antennas to infinity. 
Monte-Carlo simulations used to support our findings indicate 
that the results are useful for systems with moderate numbers 
of antennas which is consistent with other works on multi- 
antenna systems utilizing infinite random matrix theory [14|. 
Furthermore, as hardware costs decrease and transmissions 
move to higher frequencies, it will be possible to fit large 
numbers of antennas on standard mobile devices. For instance, 
at a nominal carrier frequency of 6GHz, it is possible to fit 
30 or more antenna elements on a standard laptop computer 
with wavelength separation (~ 5 cm at 6GHz). 

We start by showing that with appropriate normalization, 
the SINR of a link located at the center of a network 
converges with probability 1 to an asymptotic limit that is 
dependent on an intensity function A(r, 9) where rA(r, 9)drd9 
approximately equals the average number of points that occur 
in a small rectangle of sides dr and rd8 around the point 
(r, 9). We apply this result to power-law intensity functions 
where A(r, 9) — pr e for some e, to find an approximation 
for the SINR when N is large. This result which includes the 
uniform density as a special case (i.e. e = 0) can be used to 
study the scaling behavior of the SINR as the user density 
increases. In addition to the case where the density of users 
increases uniformly (p — > oo), this result also enables us to 
study a scaling regime where the user density "flattens" i.e. 
e -> which can be used to model networks where the user 
distribution approaches a homogeneous intensity. 

Additionally, since we make no assumptions on the distri- 
bution of the channels between pairs of antennas, our results 
can be adapted to Direct-Sequence Code-Division-Multiple- 
Access (DS/CDMA) systems with random spreading codes. 



We present the system model in Section II followed by 
the main results for general intensity functions and power 
law intensity functions in Section III. Section IV contains 
numerical results followed by a summary and conclusions. 
Outlines fof the derivations of the main results are presented 
in the appendix. 

II. System Model 

Consider a circular network of radius R centered at the 
origin. Assume that n transmitting nodes are distributed inde- 
pendently and randomly in the circle and that the i-th node is 
at location (r^, 9i) where the angle and distance are distributed 
jointly according to 

f r ,e(r,d) = rA(r,6). (1) 

Additionally, assume that n and R are related by the following 
equation: 

n= / rk{r,9)drd9. (2) 

JO Jo 

which we assume is invertible with g(n) = R and that A(r, 9) 
is such that n — > oo as R — > oo, i.e. there are an infinite 
number of nodes in an infinite radius network. 

We shall refer to these n transmitting nodes as interferers. 
Assume that a receiver is located at the origin and is com- 
municating with an additional transmitter located at a fixed 
distance rx away. We shall refer to this transmitter as the target 
transmitter. The interferers are assumed to be communicating 
with other receivers whose locations do not impact the results. 

Suppose that all nodes transmit with unit power and that the 
average power received at the origin from node-i, pt follows 
the standard inverse power-law model 

Pi = r^ a (3) 

with a > 2. Each transmitting node has a single antenna and 
the receiver at the origin has N antennas. If the N x 1 vector 
of received samples at each antenna of the receiver at a given 
time is contained in the vector y, the system can be described 
by the following linear equation 

y = HP5x + w (4) 

where the N x (n + 1) matrix H contains independent, 
identically-distributed (IID) complex random variables of unit 
variance. Note that we do not require the channel co-efficients 
to be Gaussian. P = diag(r^, a ,r^ a , • • ■ r~ a ) and x is an 
(n + 1) x 1 vector of transmit samples from the target trans- 
mitter and n interferers. w is a vector of IID noise samples 
with variance cr'jq — <j 2 Ng(n)~ a . We define the noise in this 
manner to enable an asymptotic analysis in the interference- 
limited regime as the number of receiver antennas N — > oo. 
Without this normalization, as N — > oo, the MMSE receiver 
will be able to suppress interference to levels comparable to 
the thermal noise resulting in the system no longer being 
interference-limited. Since our focus in this work is on the 
interference-limited regime, the actual value of the thermal 
noise power will not be relevant in the results. Additionally, 



we define a normalized version of the SINR j3^ which we 
show converges with probability 1 in the next section 



(3 N = r%{g(Nc))- a SlNR 



(5) 



where c = n/N . 

The receiver uses a linear-MMSE estimator to estimate the 
sample from the target transmitter, i.e. the first entry of x. The 
SINR associated with this estimator is known to equal 



SINR = hl( HPH t + (j 2 N l) h T . 



(6) 



where is the first column of H, P = diag(pi,p2, • • -p n ) 
and H is equal to H with the first column removed. 

III. Main Results 

A. General Intensity Functions 
Theorem 1: Let 

fjf" A(g{Nc)x-i, 9)l { i <x<oo} d9 



H{x) =l-x~ 



lim 

N->-oc 



f* n A(g(Nc),6)d0 



Then, in the limit as n,N,R — > oo such that c = n/N and 
(|2]l always holds, fjjf — > j3 with probability 1, where f3 is a 
unique non-negative solution to the following equation: 



-(3a 2 + l=fic 



x dH(x) 
l + xfi 



(7) 



Thus, when the number of antennas at the receiver is large, 
the SINR can be approximated by removing the scaling as 
follows: 



SINR « r T a f3g{n) a = r^ a /3g(Nc) a . 



(8) 



Note that in many cases of interest (e.g. the power law intensity 
function described below), the RHS above is does not depend 
strongly c. 

Theorem 2: If (g(n)) a is a function of lower than expo- 
nential order, i.e. for all A > 1, A~ n (g(n)) a — ► as n — > oo, 
the following hold with probability 1: 



log 2 (l + SINR) - log 2 (l + r T a (g(Nc)) a (3) -> (9) 



and 



£[log 2 (l + SINR)] - log 2 (l + r T a (g(n)) a (3) -> 0. (10) 

The latter two results are useful to characterize the spectral 
efficiency when Gaussian code-books are used. 

B. Power-Law Intensity 

Now consider an intensity function A(r, 6) — pr e for e > 
—2. This model can be used to approximate systems where 
the intensity of nodes decays with distance from the origin 
such as that depicted in Figure Q] 

Theorem 3: Under the conditions of Theorem Q] /3jy — > (3 
with probability 1, where (3 is a unique positive solution to: 



/3a 2 + 1 = —/3~ttcsc 1 



(2 + e)/3c 



2^1 1, 
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-2-e + a -2-e + la 



a 



(11) 

With e = 0, which corresponds to a uniform intensity of 
nodes, the result above is a strengthening of an earlier work 
[2 1 where it was shown that E\j3jsr] — > (3 and var{/8jv} 0. 

For a — 2 — e > 0, we can find a simple approximation 
for (3 in the expression above if we write t/n — Pnc^ = 



SINR and rj = c 2 + r - (3, we have: 



77c 2 +' a + 1 



(2 + e) 2+s /27r + 7re 
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(2 + e)r]c^ 



■ x 



1F1 1 



a-2-e 

-2-e + a -2-e + 2a 



-rjc 



a a 
Using Euler's Hypergeometric transforms yields: 

^2 , (2 + e) 2+e (2-K + Tie 

- rjc 2 +« a + 1 = -rj ° 7r esc 



a 



a 



(2 



2 + e-a / -a 

■qc 2 +' / -2-e + 2a rjc 2 +' 

■2F1 1,1, , : 

\ a 



a - 2 - e 1 _|_ ^ C 3q^ y a 1 + ^3+7 / 

Since a — 2 — e > 0, the parameters of the hypergeometric 
function above are positive implying that it is a monotonically 
increasing function of the argument. Hence, 



1 



(1 + 7iC 2 +0 



1F1 1,1, 



-2-e + 2a rjc 2 +' 



a 1 + 77CS+7 

2 + e-a 



< 



a / 2 + e 

For large c, i.e. when the ratio of nodes in the network to the 
number of antennas at the receiver is high and small a 2 : 



{2 + e) 2+e /2?T + 7re 
-rj " ix esc 



1 

a \ a 

Thus for large N, the SINR can be approximated by 

'27r + 7re^* ^N\^ 



SINR 



■ sin 



(12) 



(13) 
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and by Theorem [2] the spectral efficiency assuming Gaussian 
codebooks can be approximated as: 

27r + 7re^ 5 T 7 ^iV^j^^ 
(14) 



7 s» log 



1 



2tt 2 



■ Sill 



Re-writing the above equation in terms of N yields 



N w p4 +e (2 7 - 1) 



2+e /2tt 



■ CSC 



27T + 7T6 



(15) 



This indicates that to maintain a constant spectral efficiency 
N must increase linearly with p and r 2+£ but approximately 
exponentially with e. 
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Fig. 2. Normalized SINR vs. number of receiver antennas for power-law 
intensity with e = —0.5 and a = 3. The points are 50 random samples out 



of the simulated values of r;jv = (^U^-^ 
lines indicate the asymptotic approximation. 
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Fig. 3. Mean spectral efficiency vs. number of receiver antennas for power- 
law intensity with e = —0.5 and varying a. The markers indicate simulated 
values and the solid lines indicate the theoretical approximation of (14) . 



IV. Numerical Results 

We performed Monte-carlo simulations of network topolo- 
gies to validate the results and approximations from the previ- 
ous sections. The receiver was placed at the origin and 10000 
nodes were placed at random polar coordinates in a disk of 
radius R which satisfies © with p = 0.01. The distances were 
drawn from a density proportional to r e+1 between and R 
and the angles were uniform in (0, 2ir]. The channels between 
the antennas of a node at a given distance and the receiver were 
IID, circularly symmetric Gaussian random variables. We used 
unit transmit power and a constant noise power of 10~ 14 and 
the target transmitter was at distance tt = 10. Note that in 
this regime, the system is interference limited and the specific 
value of the noise power does not appreciably impact the 
spectral efficiency. 1000 trials of each set of parameters was 
simulated and the spectral efficiencies were computed using 
the Shannon formula. We considered a — 2.5,3.5,3,4, and 
4.5 and e = -1, -0.75, -0.5, -0.25, and 0. 

The points in Figure [2] illustrate a random sampling of 50 

points from the simulation of r/w = ff^f -^H 2+ ' r r SINR 
for each value of N, and the solid line shows the asymptotic 
prediction for a system with e = —0.5 and a ~ 3. We plotted 
a small subset of the simulated points for clarity. The concen- 
tration of points with increasing N indicates convergence of 
the normalized SINR to the asymptotic value. 

Figures |3] and |4] show the spectral efficiencies for e = —0.5 
and different values of a. Figure [3] shows the simulated mean 
spectral efficiencies (markers) and the asymptotic prediction 
of (Tl4l . In all the cases shown, the mean spectral efficiency is 
within 5% of the asymptotic value for N > 6 which illustrates 
the utility of ( fT4l for reasonably small numbers of receiver 
antennas. Figure [4] shows a random sampling of 50 trials of 
the spectral efficiency for each value of N, for a = 2.5 and 
a = 4. The concentration of points indicates convergence of 
the spectral efficiency to its asymptotic value. 

Figure [5] shows the mean spectral efficiency for a = 4 and 




10 20 30 40 50 60 70 

Number of Receiver Antennas 



Fig. 4. Simulated spectral efficiencies vs. number of receiver antennas for 
power-law intensity with e = —0.5 and different values of a. The points 
indicate a random sampling of 50 simulated values, and the solid lines indicate 
the theoretical approximation. 



different values of e. The mean spectral efficiency is within 
5% of the asymptotic value for N > 6 indicating that (fT4l 
holds for a range of different e and reasonable numbers of 
antennas. Figure |6] is a plot of N vs. e from (fT~5T > for different 
values of the spectral efficiency. Note that an approximately 
exponential increase in N is required to maintain a constant 
spectral efficiency with increasing e. Furthermore, large spec- 
tral efficiencies greatly impact the rate at which N needs to 
increase with e. 

V. Summary and Conclusions 

We have developed asymptotic technique to find the SINR 
and spectral efficiency of an interference-limited link with 
a multi-antenna receiver in a wireless network with nodes 
distributed according to a non-homogenous point process. We 
find that with appropriate normalization, the SINR and spectral 
efficiency of this link converges with probability 1 to an 
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Fig. 5. Mean spectral efficiency vs. number of antennas for different values 
of e with r T = 10, a = 4, B = 0.01. 
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(e — »■ 0) rather than increases uniformly. We find that an 
approximately exponential growth of N is needed to mantain 
constant spectral efficiency as e — >• 0. 

The asymptotic results are supported by monte-carlo simu- 
lations for which the simulated mean spectral efficiencies are 
close to the asymptotic prediction, differing by less than 5% 
when the number of antennas is greater than 6, for all the 
cases considered, confirming the the utility of these results for 
finite numbers of antennas. 

In addition to strengthening previously known results for 
homogenous networks, these results characterize the spectral 
efficiency of a link with multiple receiver antennas in networks 
with a non-homogenous distribution of nodes for which there 
are few results in the literature. 

Appendix 

A. Outline of the derivation for general intensity functions 

We accomplish the normalization of the SINR by (g(n))~ a 
by scaling the interference and noise powers by (g(n)) a . Let 
the scaled interference power due to node i be pi = (g(n)) a pi 
and the thermal noise power be a%{g{n)) a . Consider the CDF 
of pi for a given R and recall that R = g(n). 

Fp(x) = Pr(r~ a (g(n)) a < x) = 1 - F r (aT^n) 

_ n j;"°» w f r AM)y m rt 



{r<R} < 

J Q R j^rA(r,9)drd9 



jRx « r ^Q}l {r<R}drd Q 
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Fig. 6. Number of receiver antennas required to maintain constant spectral 
efficiencies vs. e. Link-length rx = 10, a = 3, B = 0.01 



asymptotic limit as the number of receiver antennas and radius 
of the network grow to infinity. 

We applied this result to networks with intensity functions 
of the form pr e for e > —2, to derive a simple approximation 
(given by (fl4] i) for the spectral efficiency. This class of 
intensity functions can be used to model receivers located 
in the center of a dense cluster of nodes and includes the 
homogenous node distribution as a special case (e = 0). 
For this model, the SINR is found to be approximately 
proportional to N/p raised to the power which indicates 
that per-link spectral efficiency can be held constant by linearly 
increasing N with density p, which has been found to hold for 
homogenous networks in 0. Comparing the e = — 1 case for 
which the SINR grows as N a and the homogenous case for 
which the SINR grows as N~z , we find that the SINR grows 
with TV at a rate significantly faster in the center of a dense 
cluster compared to a homogenous network. 

Additionally, our results enables the analysis of a new 
scaling regime where the intensity of the node process flattens 



For any positive x, as R — > oo, both the numerator and 
denominator must go to infinity and thus the limit is the ratio 
of the deriviatives w.r.t. R of the numerator and denominator. 
Observing that — > oo as R — > oo we have: 



lim Ffi(x) = 1 

R^roo F 



hC~" CrA(r,9)l {r<R} drdB 



d J R CrA(r,9)drd9 



dR 

CA(Rx- 



lim ■ 



oor 



"<i 



d6 



J^A(R,9)d9 



lim 



A{g(Nc)x-i,9)l {1<x<oo} d9 



J^A(g(Nc),9)d9 



(16) 



Let H(x) = lirrifl_j. 00 F p (x). By the main result of |fT31l 
and the Glivenko-Cantelli theorem, we have f3jy — > f3 with 
probability 1, where (3 is a unique non-negative solution to 
the following equation: 



-(3a 2 + I=f3c 



x dH(x) 
1 + xP 



(17) 



which completes the proof of Theorem Q] 



To prove the first part of Theorem [2] consider the regime 

as N — > oo: 

log 2 (l + SINR) - log 2 (l + r T a g(Nc) a p) 

= log 2 (l + r T a g(Nc) a p N ) - log 2 (l + rj a g(Nc) a P) 

= log 2 (g(Nc)- a + r T a p N ) - \og 2 (g{Nc)- a + r T «/3) -> 

where the last limit holds with probability 1 by the continuity 
of the log function. The second part can be proved using 
similar steps as Appendix E of Ifl6l . 

B. Outline of the derivation for power-law intensity functions 

The intensity function is A(r, 9) = pr e for e > — 2 which 
yields the following: 



2^ 2 



9{n) = R 



2np 



-Nc 



H(x) = 1 - l { i <x} Yirn^ 



2np 



dH{x) 



1 - 1{1<i}S 

2 + e 2+,+q 
x 



{Kx} 



dx a 
which when substituted into dT) yields: 

2 + e .2+, 



1 



/3c 



(2 + e)/3c 



o a(l+ J 9a;) 



da; — 



^{i<x}dx 



(18) 

(19) 

(20) 
(21) 

(22) 
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o + 

27T + 7T6 



1 -Ml 
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(2 + e)/3c 



Fi 1, 



a 

2 - e + a 



2a 



a — 2 — e \ ol °- 

which concludes the proof of Theorem [3] 
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